Coupling of replicate order-parameters in incommensurate multiferroics 
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The specific properties of incommensurate multiferroic phases resulting from the coupling of order- 
parameter rephcates are worked out using the illustrative example of iron vanadate. The dephasing 
between the order-parameter copies induces an additional broken symmetry phase corresponding to 
the lowest symmetry of the system and varies critically at the transition to the multiferroic phase. 
It reflects the temperature dependence of the angle between paired spins in the antiferromagnetic 
spiral structure. Expressing the transition order-parameters in terms of spin-density waves allows 
showing that isotropic exchange interactions contribute to the stabilization of the ferroelectric phase. 

PACS numbers: 77.80.-e, 61.50.Ah, 75.80.-|-q 



One of the distinctive features of magnetic multiferroic 
materials is the existence of an incommensurate ferro- 
electric phase induced by the coupling of two antiferro- 
magnetic order-parameters. In multiferroic compounds, 
such as TbMnOa,^ Ni3V20s^ or MnW04^ the two order- 
parameters display different symmetries, i.e. they trans- 
form as distinct irreducible representations (IR's) of 
the paramagnetic space-group i^i^ A different situation 
has been pointed out in the pyroxene NaFeSi206,^ in (772 



which the onset of the ferroelectric phase results from 
the coupling of two dephased copies of the same order- 
parameter. The possibility of inducing an additional 
broken-symmetry phase by replication of a single order- 
parameter was shown^ to reflect the existence of an effec- 
tive continuous symmetry which is broken at the transi- 
tion to the ferroelectric phase, expressing the continuous 
rotation of the phason mode associated with the incom- 
mensurate spin- wave order-parameter. 

Phase transitions induced by the coupling of replicate 
order-parameters display a number of specific features, 
such as the stabilization of additional broken-symmetry 
phases and a critical dependence of the dephasing be- 
tween the order-parametersj^ These properties are illus- 
trated here by the theoretical description of the tran- 
sition to the multiferroic phase observed in iron vana- 
date FeV04 which, due to the low triclinic symmetry of 
the paramagnetic phasCf^-iS constitutes an unambiguous 
proof of replication mechanism. Expressing the transi- 
tion order-parameter in terms of magnetic spin waves 
allows showing that the order-parameter dephasing cor- 
responds to the angle between spins forming pairs in the 
antiferromagnetic spiral structure, and that isotropic ex- 
change interactions contribute to the emergence of the 
polarization in the multiferroic phase. 

Starting from the paramagnetic symmetry Pll' of 
FeV04, a single 2-dimensional IR (T) can be constructed 
in the general direction 4 = (0.222, 0.089, 0.012) of the 
triclinic Brillouin-zone reported experimentally the ma- 
trices of which are given in Table HI A symmetry anal- 
ysis of the corresponding 2-component order-parameter 
{r]i = pie^^^ ,rjl = pie~^^^) shows that a single second- 
order transition to an incommensurate non-polar phase 



of point-group 11' (phase /) can be induced by F, which 
coincides with the antiferromagnetic phase observed be- 
tween Tp = 22K arid Tp = IbK in FeV04.^ Since the 
propagation vector k is nearly temperature- independent, 
a further symmetry breaking mechanism to the ferro- 
electric phase stable below Tp cannot be obtained if 
not assuming that a coupling takes place at Tp be- 
tween the order-parameter (?7i,7?x) and a dephased copy 



P2e''\V*2 



P2e~ ) having the same F-symmetry, 
as already noted in ref. [8] The free-energy expressing 
the coupling between the two copies can be written: 
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where 6 — 9i — 62- Fig. [T] shows the phase diagram de- 
duced from the minimization of F with respect to pi , p2 
and 6. Below the range of stability of phase I [9 — mr), 
the phase diagram contains an additional phase // of 
magnetic point group 11', corresponding to 9 ^ mr, 
which displays a spontaneous polarization P in general 
direction, as observed below Tp in FeV04^'^ The topol- 
ogy of the phase diagram of Fig. 1 is specific to replicate 
order-parameters, the symmetries of phases / (11') and 
// (11') representing, respectively, the minimal symme- 
try group induced by F, defined as the invariance group of 
the general direction in the representation space, and the 
lowest symmetry group induced by F -t- F corresponding 
to the kernel of the homomorphism of the paramagnetic 
space group Gp on F. 

The straight thermodynamic path shown in Fig.l 
can be parametrized as ai — x{T — Tp) + 4^ and 
bi = ^(T — Tp) -\- w, where the constant parameters 
X, C ^iid w characterizing the path fulfil the condition 
xC{Tp-Tp)^+4.^{Tp-Tp)+3cl ^ 0. It yields the prop- 
erty that the order-parameter amplitudes pi — (~^)^^^ 
and p2 = (~^)^^^ vary non- critically at Tp, while the 
dephasing 9 plays the role of the critical transition order- 
parameter, varying below Tp as: 

.^n0=(g + ^)'''^/7W (2) 
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FIG. 1: (Color online) Phase diagram associated with the 
free-energy F given by Eg. All curves correspond to 

second-order phase transitions. The arrow represent the ther- 
modynamic path used for calculating the dephasing 6 given 
by Eq. ©. 



with opposed polarizations corresponding to opposed de- 
phasings ±9. 

Application of a magnetic field does not give rise 
to spectacular magnetoelectric effects in FeV04, since 
only trivial couplings between the order-parameter and 
the magnetic field are allowed by the low symmetry of 
the system. Thus, biquadratic couplings of the order- 
parameter to the magnetic free-energy fiijM^M^ yield a 
renormalization of the ai,bi, and Ci coefficients in Eq. 
(H)) explaining the observed shift of Tp to lower temper- 
atures^ with increasing magnetic field. Standard mag- 
netoelectric couplings P'^M^ are responsible of the 
progressive reduction and suppression of the multiferroic 
phase reported at high fieldi^ A different manifestation 
of magnetoelectricity can be found in the toroidic ef- 
fects which can be predicted in the ferroelectric phase 
of FeV04. Under applied field B the toroidal free-energy 
is: 



(4) 



where T' = k are the induced toroidal moment com- 



TABLE I: Generators of the irreducible representation F of ponents and 



-\piP2sin6{g ^Y^aij the magneto- 



the Pll' space group for k = {k^,ky,kz). I is the space toroidal susceptibility components. It yields: 
inversion, T is the time reversal and R = [na + pfc + qc). 
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This property differs from the standard behaviour found 
in multiferroic phases induced by the coupling of order- 
parameters having different symmetries where the de- 
phasing = 7r/2 is constant and the order-parameter 
amplitudes pi vary critically with temperature. We 
show below that the critical variation of 9 can be in- 
terpreted as the temperature dependence of the angle 
between paired spins which are parallel in phase I and 
become non-parallel in phase //. The equilibrium po- 
larization deduced from the coupling free-energy F = 
aiP^ pip2sin9 -\- XijP^P^ , gives in triclinic basis: 



P' = -^iX-^y'ajPip2Sin9 



(3) 



Close to Tp 

1. ^T/Td -L /I ^2 ^ 



one gets 

1/2 



AWTf ~ T with A* = 
{X^^Y^aj. This power law and the 

correlated Curie- Weiss like sharp peak of the dielectric 
permittivity e are consistent with the dependence on tem- 
perature of P{T) and e{T) observed on single crystal^ 
and polycrystallineS. FeV04 samples. It confirms the hy- 
brid pseudo-proper character^ of the ferroelectric tran- 
sition which behaves critically as a proper ferroelectric, 
although the small value of | P | found below Tp (6- 
12 nCm~^) is of the order of magnitude of an improper 
ferroelectric. Since inversion symmetry is broken at the 
transition the system displays two ferroelectric domainsil 



r' = K^^Tp - T 



(5) 



with K 
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-5(x«^ + 4^j' {g-'T^-jBK Eq. ® 
shows that at constant field the toroidal moment varies 
critically as the spontaneous polarization at zero fields. 
Since all components of the magnetotoroidal tensor are 
non-zero the induced toroidal moment should not orient 
along B. 

To gain insight into the microscopic spin states form- 
ing the magnetic structures, one can express the transi- 
tion order-parameter components in terms of spin-waves. 
The triclinic centred unit cell (Fig. [2]) in which the inver- 
sion centre is assumed to coincide with the centre of the 
cell, contains 6 magnetic atoms in positions (1,2,3) and 
(T, 2, 3). Denoting j the spin of the i atom in the cell 

centred at Rnpq — na pb -\- qc, with m=lorr7i = — 1 
for atoms (1, 2, 3) or (1, 2, 3), respectively, one can write: 



s™ - = - 



[7, -I- imWi 



U* - imW* 



„ik-R„ 



-ik-R„ 



(6) 



where Uj and Wj are spin-waves transforming as V for 
each j and each space component. It yields 18 spin- waves 
which are copies of 77 and 77*. Eq. ^ can be written 

^/c ■ Rnpq ' 



npq,i 



in real form as: 
mp'^ ^sin {k ■ Rnpq ^ 
1,2,3 with C/f = p^^e*' 



j a^, where = a,b,c for /i = 
md — P2 ^e'^^.i ^ In phase 

/ all Uj and Wj are in-phase and one can chose a domain 
where they are all real. In this domain the inversion 
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FIG. 2; Positions of the Fe atoms in the trichnic unit-cell of 
FeV04. Atoms forming the pairs (1,1), (2, 2) and (3,3) are 
related by inversion. 

centre is located at the centre of the unit-cell {n = p = 
q ^ 0) and the spins take the simpler form: 

^pq.i = UiCOs{k ■ Rnpq) " mWiSin{k ■ Rnpq) (7) 

analogue to Eq. ^ of ref. [7] In order to estimate the 
dephasing 6 between the order-parameter replicates let 
us consider the six spins contained in the unit-cell with 
n — p = q = 0, namely with j = 1, 2, 3. In the non- 
polar phase / the inversion centre present at the centre 
of the cell imposes that = s~ for all j, although sf ^ 
S2 ■ In the polar phase // the loss of inversion 

centre yields non-parallel spins for each pair, and one 
has: 

_ ^ 2pi; jsine!^ . (8) 

where the angles 62 j are replica of 9, all sinOi^ j varying 
critically with temperature at the I ^ II transition as 
sine (Eq. dU). 

It has been suggested in ref. (7) that Dzialoshinskii- 
Moriya (DM) interaction may not be the sole active 
mechanism at the transition to the multiferroic phase 
of FeV04. This property can be verified by considering 
only one of the two-atoms sublattices, as for example the 
sublattice involving the (1, 1) pair, and neglect its inter- 
actions with the two other sublattices. Since the atoms 
of a single pair are related by inversion, no DM Inter- 
action can take place. We will now show that a polar- 
ization can be induced in the sublattice when assuming 
only exchange interactions between the two spins. This 
will demonstrate in turn that purely exchange contribu- 
tions (one for each sublattice) participate to the total 
polarization in FeV04 in addition to the DM and other 
frustration contributions. 



Denoting s^'^ the /i-component of the spin = 

s^pqim = ±1) in the unit-cell n = na + pb + qc the 
Hamiltonian of the system becomes: 

^ = E E h^:!^r''^''s^^'4''' (9) 

1 ira,m' ,11,11' j —m,—m',n,n' j i ■ 

wiiere /i- — due to mversion symme- 

try. At low temperature the spin modulus is frozen and 
only the orientations of the spins can vary. Neglecting 
quantum effects and higher harmonics the ground state 
of the system is obtained by minimizing H restricted to 
the first harmonic with respect to the spin directions. 
Inserting ([7]) into (O one has for ^ 0- 

H= {U" + imW''){U'''* -im'W'''*)A';^';^, +CC 

(10) 

where A^^, = E„- /»:^"'.,.,;e'^^(«"-^") = ^'l^^^. 
When only isotropic exchange interactions are involved 

"n-n' ~ ''n-il'" ^^^^ ^m,m' — ^m,m"J 

Therefore H reduces to: 

H = aU - U* + I3W - W* +j{U -W* + W -U*) (11) 

where a = 2{Aii + Ai^i + A^n + A-1-1), l3 = 2(Aii - 
Ai-i - A-u + A_i_i) and 7 = 2i{-Ai^i + A^n). 
Writing the complex vectors as U — ui + iu2 and 
W — wi + iw2, the conditions of constant spin modu- 
les II 11= So read: 

u\ + ~ + w\ ~ So/2, ui • zZ;2 + "1*2 • wi = 

Ui ■ U2 = Wi ■ W2tUi ■ Wi — U2 ■ W2 (12) 

and H becomes: 

H = a{ul + ul) + I3{wl + wl) + 2-i{uiWi + U2W2) (13) 

Phase / corresponds to ui = ±U2 and wi -f ±W2- 
Accordingly, the exchange Hamiltonian has the con- 
stant value Hq = ^{a + p). The stability of this 
phase can be tested by varying slightly the four vec- 
tors ui,U2,wi,W2- Expanding the constraints given by 
Eqs. ([T2I) to the first-order shows that only five indepen- 
dent parameters {uj'^,Suix,y,z,3u2x.y,Swiy,Sw2y), where 
= Sulx + Su2x + 2(5^^^, determine the stability of 
phase /. At the second order the Hamiltonian reads: 

H = Hq + a{u? + Sujy + 5uly) + I3{uj^ + Swjy + 

S^2y) + 2^[6uiySu2y + 6wiy5w2y] (14) 

The transition from phase / to the multiferroic phase // 
associated with the wave vector kc occurs when = 
and I a{kc) \ = \ 7 | or | (3{kc) \ — \ 7 |. These equations 
give, respectively, the equilibrium value of kc in phase 
/ and the stability limits of phase //. In order to show 
that these limits exist and that a polar phase can actually 
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be stabilized below phase /, we consider the situation 
where only interactions between first-neighbouring cells 
are involved. In this case the non-zero coefficients in Eq. 
([TUl) are h]'^, /i^^ , and one has: 

a{kc) = (2/i" 4- hl-^ + h^^^)coska ■ a 
+ {a -^b) + {a^c) 

= {2hl^ - hi-' - h^'')cosk, ■ a 
+ {a -^b) + {a^c) 

7(4) = 4(/ii-i - /i- ")sinfc, • a 

+ {a^b) + {a^c) (15) 

If hi' > and hl'^ < Hq is minimum at fc^ = (^, 0, 0) 
and a becomes hl^' - {2h}} + hlr' + /i-") + {2hl' + 
hl~' + /ig") + {2hy + hi-' + hz''), with an analogous 
expression for (3 and 7. When the h'~'^ and h^ coeffi- 

cients vary, a ± 7 and /3 ± 7 take positive and negative 
values. Thus the non-polar phase I and polar phase II 
can both be stabilized by isotropic exchange interactions. 
Considering interactions between spins belonging to 
more distant unit-cells would stabilize commensurate or 
incommensurate, polar or non-polar phases depending 
on the values of the /i™'™, coefficients of the 

exchange Hamiltonian. However, since no assumption is 
made on the actual values of the interaction coefficients, 
one cannot exclude that in h^'-i^' space, the region of 

a,h,c 

the phase diagram in which the polar phase is stable 
would actually be forbidden. 

In summary, the specific features of multiferroic phases 
induced by the coupling of replicate order-parameters 
have been illustrated by the phenomenological descrip- 
tion of the sequence of triclinic phases observed in 
FeV04. The replication mechanism yields the onset 
of an additional lower symmetry ferroelectric phase 



in the phase diagram, the dephasing between the 
order-parameter replicates varying critically at the 
ferroelectric transition. Considering a single sublattice 
of FeV04 the stabilization of the ferroelectric phase has 
been shown to result from purely isotropic exchange 
interactions. It would be of interest to verify if this 
property is specific to mutiferroic systems involving a 
replication mechanism. 

The possibility of stabilizing additional phases by cou- 
pling two dephased order-parameters having the same 
symmetry, i.e., transforming as the same irreducible 
representation of the parent phase, was first noted 
in superconductors^^ which involve a breaking of the 
continuous gauge symmetry. Another example of repli- 
cation mechanism was described in ferroelectric liquid 
crystals^"* in which the continuous rotational symmetry 
is broken at the transition to the ferroelectric phase. 
Although structural and magnetic phase transitions 
correspond to discrete symmetry breaking processes 
which do not give rise to additional phases by order- 
parameter replication, transitions to incommensurate 
structural or magnetic structures permit a realization 
of the replication mechanism, as shown in the present 
work for the multiferroic transition in FeV04. This is 
because the thermodynamic features of the transition 
are driven by an effective continuous group generated 
by the images of the crystallographic translations in the 
order-parameter space: The discrete translation group 
broken at the transition to the incommensurate structure 
has a dense image in the order-parameter space acting 
as the continuous rotation group of the phason. In other 
words, for incommensurate wave- vectors the infinite 
discrete group of translations has the same effect on the 
form of the transition free-energy as a continuous group. 

The authors are grateful to Laurent Chapon and 
Dmitry Khalyavin for helpful discussions 
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